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Abstract 

This paper studies two related questions in high energy onium-onium scattering: the 
probability of producing an unusually large number of particles in a collision, where it is 
found that the cross section for producing a central multiplicity proportional to k should 
decrease exponentially in \fk. Secondly, the nature of gluon (dipole) evolution when dipole 
densities become so high that saturation effects due to dipole-dipole interactions become 
important: measures of saturation are developed to help understand when saturation 
becomes important, and further information is obtained by exploiting changes of frame, 
which interchange unitarity and saturation corrections. 



*This work is supported in part by the US Department of Energy under Grant DE-FG02-94ER40819 
^This work is supported by the UK Particle Physics and Astronomy Research Council 



1 Introduction 



The Balitsky, Fadin, Kuraev and Lipatov (BFKL) equation |T], ||, ||] is meant to describe 
the energy dependence of certain high energy hard processes. For the BFKL equation to 
be applicable certain conditions must be met: (i) The process, or part of a process, to be 
described should have only a single hard transverse momentum scale. Processes having more 
than one hard scale may have significant contribution from the Dokshitzer, Gribov, Lipatov, 
Altarelli and Parisi (DGLAP) equation [||, ||, |(| which governs scale dependence in QCD. 
Certain observables (e.g. f7|) are known to satisfy this criterion though they have not yet 
been well measured, (ii) The rapidity of the hard process, related to the centre of mass 
energy E and the hard transverse momentum scale k by Y = In E 2 /k 2 , should be large. 
The BFKL equation includes all terms (a s Y) n in the coupling but neglects factors of a s not 
accompanied by factors of Y. There is a vigorous effort [||] to systematically include the first 
non-leading terms, an a s unaccompanied by a Y and this should improve the accuracy and 
reliability of the BFKL equation. 

The BFKL equation should describe high energy processes having a single hard scale in a 
limited range of rapidity. When Y becomes too large, unitarity corrections become important 
and slow the rate of growth with energy. At even higher values of energy one expects parton 
densities to become so large that the whole perturbative QCD picture breaks down leading to 
a new regime of weak coupling but high density non-perturbative QCD. The goal in small-x 
physics is to understand, experimentally and theoretically, the onset of BFKL behaviour, its 
modification due to unitarity corrections as one proceeds to higher energies and finally the 
nature of the transition to the high density (saturation) regime at even higher energies. 

From a theoretical point of view the simplest process to study is high energy heavy 
quarkonium-heavy quarkonium scattering. If the onium is sufficiently heavy the radius of 
the onium is small enough to set the single hard scale of the process. Any understanding of 
BFKL dynamics and beyond gained in studying heavy onium collisions should be applicable 
to the more experimentally accessible single hard scale processes mentioned above. In the 
large N c limit of QCD a dipole picture of high energy heavy onium-heavy onium scattering 



has been developed which is equivalent to the usual BFKL description |j, 10, [Ll|, [T^, [L3|]. In 
this dipole picture the light-cone wave function of a heavy onium is viewed as a collection of 
colour dipoles. The dipoles are made up of the quark part of a gluon (or the heavy quark 
itself) and the antiquark part of a different gluon (or the heavy antiquark). High energy 
elastic scattering in the BFKL approximation proceeds by two-gluon exchange between a 
right-moving dipole in the right-moving onium and a left-moving dipole in the left-moving 
onium. 



1 



Perhaps the main advantage of the dipole picture is that the equation governing the 
number and distribution of dipoles in the light-cone wave function of the onium is a branching 



equation which allows Monte Carlo methods to be used. In a recent paper [14], one of us has 
studied high energy onium-onium scattering from energies where the BFKL approximation 
is valid into the higher energy regime where unitarity limits, for a fixed impact parameter of 
the scattering, are reached. In the centre of mass system unitarity corrections are simply the 
independent scattering of two or more dipoles of the right-moving onium with two or more 
dipoles of the left-moving onium. Because of the large number of (small) dipoles in an onium 
wave function unitarity corrections become important when the dipoles are relatively dilute 
in the onium so that independent multiple scatterings are the natural leading corrections to 
the BFKL approximation. 

The purpose of the present paper is to study two related questions, the probability of 
producing an unusually large number of particles in an onium-onium collision and the nature 
of gluon (dipole) evolution when dipole densities become so dense that saturation effects 
due to dipole-dipole interactions become important. At a given energy dipole densities in 
an onium wave function are largest not in a typical configuration but rather in rare high 
multiplicity fluctuations. Thus attempts to measure saturation will naturally focus on the 
high multiplicity tail of particle production. 

One cannot directly study particle production in the dipole formalism. In the dipole 
picture the t = light-cone wave function is constructed. This is sufficient to describe an 
elastic scattering, but it is not sufficient to describe an inelastic event where one must follow 
the time evolution of the wave function to a large positive time when particles (or jets) are 
produced. However, using the Abramovskii, Gribov, Kancheli (AGK) cutting rules jljj one 
can relate particle production and forward elastic scattering involving one or more dipole 
scatterings. While there has been no detailed proof of the AGK rules in QCD it is likely that 
they are correct and we assume that this is the case in the present paper. 

We calculate the cross section for k cut pomerons when k takes values far above average. 
The multiplicity of produced particles, in the central unit of rapidity, should be proportional 
to k. We begin by studying a toy model which obeys dipole evolution but where the absence of 
transverse dimensions makes analytic solutions to the multiplicity distribution possible. Our 
main result for the toy model is given in eq. (^) where is seen to decrease exponentially 
in \fk at large k, a tail much longer than that given by KNO scaling We have not 

been able to solve the large k behaviour for the k cut pomeron contribution to a zero impact 
parameter collision in QCD, do"fc/d 2 r(r ~ 0). However, Monte Carlo results obtained using 
OEDIPUS (Onium Evolution, Dipole Interaction and Perturbative Unitarisation Software) 
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fll7f indicate a very similar /c-dependence to that found in the toy model with an e - A Vk+n 
dependence at large k. The same dependence, but with different values of A and n, also well 
describes onium-onium collisions integrated over impact parameter. The fact that the high 
multiplicity tail extends so far encourages one to believe that high parton density systems 
can be produced at reasonable accelerator energies if the right signal is found to trigger on 
the interesting physics associated with saturation. 

In section 3, we turn to a study of saturation Jig], beginning with defining measures of 
saturation. Saturation effects show up as a modification of the simple branching picture of 
dipoles in constructing the square of the onium wave function. Thus, a non-independence of 
the dipoles in an onium wave function is a measure of saturation. However, it is not known how 
to estimate the interaction of dipoles in an onium wave function directly. One expects that 
the interaction of a slightly right-moving dipole with higher rapidity right-moving dipoles in 
a right-moving onium should be similar to that of a slightly left-moving dipole with the right- 
moving onium. This latter interaction is just the scattering amplitude for a simple dipole 
with an onium state and is straightforward to calculate. Thus one measure of saturation 
of a light-cone onium wave function is the interaction probability with an opposite moving 
elementary dipole. This measure of saturation is however time consuming to evaluate. A 
more efficient measure is the overlap that a given dipole in the onium wave function has with 
all the other dipoles in the wave function. Happily, it turns out that these two measures are 
strongly correlated so that the simpler overlap measure can be used. 

The toy model discussed earlier also furnishes a useful model of saturation. Saturation 
is expected to slow the growth of the number of dipoles. By comparing the cross section 
calculated in frames where the two incoming particles share the incoming energy differently, 
a consistency condition, requiring that the cross section be frame independent, gives sufficient 
information to determine dipole branching as a function of rapidity. While at low energy the 
rate of parton evolution grows exponentially in rapidity, as in BFKL evolution, that rate 
reaches a constant (saturates) at high rapidity as shown in eq. (16). 

The toy model can be used to study the intriguing, and still not well understood, rela- 
tionship between unitarity limits and saturation. In QCD the centre of mass onium-onium 
scattering, at a fixed impact parameter, reaches the unitarity limit at energies well below that 
at which saturation effects begin to appear. However, if onium-onium scattering is viewed in 
a frame where, say, the left-moving onium has a small rapidity, the scattering looks exactly 
like the measure we have earlier discussed for saturation since now the left-moving onium is 
simply a single dipole. Thus the multiple scatterings in the centre of mass system must show 
up as saturation effects in a system where one of the onia carries almost all the energy. This 
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can be most sharply seen in the toy model by considering a toy onium scattering on a toy 
nucleus consisting of N onia. Here it is found that saturation effects found earlier are exactly 
what is necessary to make onium-nucleus scattering the same in either the frame where the 
onium is at rest or in the frame where the nucleus is at rest. 

In QCD it is not possible to do such complete analytical calculations, however, we have 
found a phenomenological modification of dipole branching eq. (33) which makes onium- 
onium and onium-nucleus scattering reasonably frame independent. One particularly sharp 
way to see the necessity of saturation effects and their frame dependence is to compare the 
S-matrix, at zero impact parameter, for onium-onium scattering in the centre of mass frame 
with that obtained in a laboratory frame at very high rapidities. At rapidities large enough so 
that the S-matrix has become very small the evaluation of the S-matrix neglecting saturation 
is determined by the very low multiplicity fluctuations of the wave function. We find, on 
theoretical grounds and from numerical simulations, that S ~ e - c (^- y o) 2 However, without 
saturation effects c is a factor of 2 larger in the lab frame than in the centre of mass. Saturation 
effects will slow the growth of dipoles in the laboratory frame so that the typical wave function 
configuration agrees with the centre of mass calculation coming from rare fluctuations. This 
is exactly what happens in the toy model as can be seen by comparing eqs. (^) and (p8|). 

At large impact parameters, r, we find that the property of the centre of mass frame, that 
saturation corrections set in much later than unitarity effects, starts to break down. This 
is brought about by several factors. At large r, the amplitude has a significant contribution 
from asymmetric onium-onium configurations, for example where one onium has a only central 
dense cluster of small dipoles (of order of the onium size): for there to be an interaction the 
second onium must have small dipoles at the relevant large impact parameter, where they are 
likely to be much more dilute. This kind of scattering looks similar to a lab frame collision, 
and hence the saturation and unitarity corrections will be similar. A second point is that one 
might expect a distribution of large dipoles which looks relatively dilute, to have only small 
saturation corrections. However, the sequence of evolution to produce this distribution will 
often have passed through a stage with a dense collection of small dipoles, whose evolution 
would be significantly altered by saturation effects, reducing the production of large dipoles 
in the first place. Overall, in the centre of mass frame, at large impact parameters, these 
effects combine to give saturation corrections of the same order as the unitarity corrections. 
Fortunately, the total scattering cross section is not too strongly affected by these large impact 
parameters, and it is still reasonable to neglect saturation effects in its calculation. 
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2 Cuts 



2.1 The AGK cutting rules 



It is most useful to start from the AGK cutting rules 15] as expressed in [19|. The cross 



section, involving the exchange of exactly k cut pomerons can be expressed as follows: 

oo 

vk = Y,(- 1 Y~ krc k F(u) > C 1 ) 

v=k 

where is the amplitude for the exchange of v pomerons. In this notation, the total cross 
section at is 

oo 

G t = lY J (rV) v ~ X F {y) - (2) 



One knows from previous work [11, 14] that the terms of the sum in eq. (^) diverge as a 
factorial, and therefore so will the terms in eq. (|l]). The reason for this divergence is that rare 
configurations of the onia contribute increasingly large amounts to the amplitudes for higher 
pomeron exchange. The solution is to sum multiple pomeron amplitudes before averaging over 
onium configurations. Applying this method here one obtains the differential cross section for 
obtaining k cut pomerons. The v pomeron amplitude for a fixed impact parameter r between 
the two onia (of sizes b and b', evolved to rapidities y and Y — y) is: 

FM(r, b, b', Y, y) = 1 ]T P 7 (r , b, y)Py (r + r, b', Y — y)(/ 7 ,y (3) 
• 7> y 

where / 7 y is the two gluon exchange interaction amplitude between a pair of configurations 7 
and 7' which occur with probabilities Py(ro, b, y) and Py (ro + r, b', Y — y) from onia situated 
at ro and ro + r. Substituting this into eq. ([[]) gives 

^ = ^E P 7(ro,b,y)Py(ro + r,b / ,y-y)(2/ 7i y) fc exp(-2/ 7i y). (4) 
• 7) y 

The cross section for obtaining no cuts is the difference between the unitarised cross section 
and Y^k=i a ki which comes out as: 



^ = E^(ro, b, y)Py (r + r, b', Y - y) (l - e"A.y ) 2 . (5) 
7,7' 

In the work that follows, y will be set to Y/2, so that we will be examining the number of 
cuts at a central rapidity. 
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2.2 Toy model calculation 



As was the case for the unitarised amplitude [11], the evaluation of eq. (0) is very difficult 



to do analytically, because it requires a detailed understanding of the distribution of dipole 
configurations for an onium. However some of the features of multiple cut cross sections can 



be extracted by examining a "toy" model which has no transverse dimensions [11]. 

The main feature of the toy model is that the probability of obtaining a large number n of 
dipoles has an exponential distribution, P n ~ e~ n / fJ, /[i, where \x is the mean number of dipoles 
(proportional to e^ ar ~^ y for an onium evolved to a rapidity y). One can show that in the 
case with transverse dimensions, the probability of obtaining a large density of dipoles also 
has a distribution which is an exponential in the density [14]. It is this similarity which leads 
to the correspondence between the toy model and the case with two transverse dimensions. 
In the toy model, the cross section for an interaction with k cut pomerons is 

ak = \\ E( 2 / a l^) /c exp(-2/a|nm)P m P ri , (6) 

m,n 

where the interaction amplitude between two dipoles is defined to be /cr| (equivalent to /~y 
above). There are two limits which can be taken. Let F^ 1 ' be the amplitude for exchange of 
one pomeron. If kF^ = kfagfi 2 <C 1, then unitarisation effects (the exponential factor) are 
unimportant, and the evaluation of the cuts cross sections is very similar to the evaluation of 



the multiple pomeron amplitudes 11]: 



a k ^kl(2FW) k e- 2k2FW . (7) 

The opposite limit is when unitarisation effects are important, fcir(i) = kfa 2 s n 2 > 1 (and 
also k/F^ 3> 1). Approximating the sums by integrals, and changing variables to n = ^/ue^ 
and m = \/ue~^ gives 

Ok - 477(2/4) fc / dud^u k e- 2fa2 s u e- 2{cosh ^^. (8) 
[i z k\ J 

Performing the £ integration using the saddle point approximation, it is clear that the integral 
is dominated by a region where m and n are similar. Performing the second integration gives 

1 (2FW* 3 \ 1/ * , . _.. , 

a k ^ TTFTTT L ex P ~\ ■ ( 9 ) 




2F(!) \ k J y I V 

The unusual feature of this result is the square root dependence on k in the exponential. 
Given that the number of cuts is related to the final state multiplicity, one might have ex- 
pected the usual KNO type [jl^] exponential distribution. The origin of the yk behaviour 
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is that, in the region of configuration space contributing to the fc-cut amplitude, the fac- 
tor (2/a|.nm) fc exp(— 2fat 2 s nm) is the most rapidly varying and fixes the dominant values 
of n, m oc \fk. The probabilities associated with the relevant configurations just fix the 
normalisation, with P{\fk) ~ e~^. 

2.3 Monte Carlo results for differential cut cross sections 

One is most likely to find a similarity between the toy model results and the realistic situation, 
by examining the differential cut cross sections, which eliminates two of the transverse degrees 
of freedom. In [14] correspondence between the toy model and Monte Carlo results for zero 



impact parameter has already been established for the multiple pomeron exchange amplitudes. 
The main difference was that where the quantity v occurred in an expression relating to v- 
pomeron exchange, it had to be replaced with {v + n) where n was a constant which originated 
from the details of the dynamics in the transverse dimensions. One might expect a similar 
effect when looking at cut cross sections, and one can examine the differential cross section 
for small impact parameters for behaviour of the form 

^*(r ~ 0) ~ exp[-A{k + n y + B] (10) 
d^r 

where A, B, n and 7 are constants to be determined. The first stage is to ensure that one 
really does have 7 = 1/2. By defining = — lnd 2 ci/ c /d 2 r and plotting 

^4-1 — ^-1 2(k + n) „. ,,, 
r fc = — k -j— = -\ 1 + 0(l/(k + n)) (11) 

(essentially the ratio of first and second derivatives) one can eliminate A and B. Figure [l] 
shows this quantity as a function of k, determined using the OEDIPUS Monte Carlo simulation 
[17]. For k sufficiently large, the linear behaviour is very clear, and the slope corresponds 



to 7 = 1/2. The value of the intercept with the k axis yields n, which has a value of 10. 
This is somewhat larger than the corresponding value found for multiple pomeron exchange 
(y — > v + n, with n ~ 3). The difference can be explained by arguing that the multiple 
pomeron and multiple cut cross sections are sensitive to the dipole distributions (and the 
manner in which they differ from an exponential) in different ways. 

The other property to examine is the variation of A with Y . One clearly cannot expect 



A to be exactly w2/i ? ( 1 )(0), however it might be reasonable for it to be proportional to 



2/F(!)(0). Tabled gives A for two values of Y, indicating that to within the accuracy of 
the Monte Carlo results, this proportionality does hold. 
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Figure 1: The upper plot shows for fixed impact parameter: a test that the cut cross 
sections show the behaviour corresponding to eq. ( p~Q| ) (actually averaged over a range of 
r = — ► 0.66); the lower plot shows the multiple cut cross sections together with a fit of the 
form eq. ( JlOj ) (for the range r = 0.2b — > 0.36). 



8 



Y 


n 


A 


F«(0) 


A/y/2/FW{0) 


10 


7 


4.45 


0.436 


2.08 


14 


10 


2.21 


1.62 


1.99 



Table 1: Parameterisations of cut distributions at central impact parameters. 



The dependence of the distributions on n is relatively weak which is why it is not determined 
more accurately. Table |] therefore confirms that for the zero impact parameter cross section 
the distribution scales as suggested by the toy model. 

2.4 Monte Carlo results for the integrated cut cross sections 

One can argue that events with k cut pomerons will have a multiplicity proportional to k. 
In an experimental situation, one might be able to determine the cross sections for various 
ranges of multiplicities, but this will correspond to cut cross sections integrated over impact 
parameter, for which there will not necessarily be any correspondence with the toy model. 
However figure |2| shows that a fit of the form eq. (10) 



a k ~exp[-A(k + ny + B] (12) 



works remarkably well. The value of n turns out to be considerably smaller than for eq. (|10[): 
this is associated with the fact that at larger values of r the value of n decreases. 

It is not clear why the exp(— Ay/k + n) behaviour also applies to the integrated cross 
section. The results do however suggest that the variation of the coefficient A with the 
rapidity is not as simple as in the toy model (or in the fixed impact parameter case). The 
numerical results for the variation with rapidity of the fraction of the cross section coming 
from different numbers of cuts is shown in figure ||. 

2.5 Spatial distribution of cut cross sections 

Just as high numbers of pomeron exchange tend to be dominated by small impact parameters, 
one can expect events with large numbers of cuts also to be relatively central, because they are 
dominated by dense configurations, which arise most commonly at small impact parameters. 
This is illustrated in figure |||, which shows a^(r > R)/a k , the fraction of the k-cut cross 
section which comes from impact parameters larger than R, for various numbers of cuts, as a 
function of R. 

One sees for example that about half the one-cut cross section is coming from a region 
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Figure 2: The upper plot shows for the integrated cross sections; the lower plot shows the 
multiple cut cross sections together with a fit of the form eq. ([To|). 
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Figure 3: The fraction of the total cross section coming from different numbers of cuts (indi- 
cated by the labels next to the curves) as a function of rapidity. 




Figure 4: The fraction of the multiple cut cross section which comes from impact parameters 
larger than R for a variety of numbers of cuts. At this rapidity, the mean number of cuts is 
~ 1.3 
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r < 3.26, whereas the first half of the 5-cut cross section comes from r < 1.46. There is of 
course a limit to how central the many cut cross section can become, since any cross section 
is bound to be spread out over a region at least the size of the onia. The no cut cross section 
is also quite central: like the many cut cross sections, it is particularly sensitive to dense 
configurations. Note that all these cross sections have significant tails at large r. This is 
to be contrasted with the 5-pomeron curve, where the tail dies off quickly: the difference 
arises because of the unitarity corrections that are associated with multi-cut cross sections 
but not multi-pomeron contributions, and which limit the relative contributions from very 
dense configurations. 



3 Saturation 



In the calculations performed to determine the unitarity corrections [11, 14], it has been 
assumed that saturation of the wave function can be neglected. The justification for doing this 
is that in a situation where both wave functions are evolved to Y/2, unitarisation corrections 
will set in when 1, while saturation corrections for a particular dipole will be 

of the order a|e( ap ~ 1 ' 5 '' 2 . Two problems can arise with this argument. The first relates to 
the question of whether as is small enough. In the calculations presented here, as ~ 0.2, so 
when unitarisation is important one might expect 0.2. Since the conditions 

for saturation can only be guessed to within a factor of order 1, one is potentially close to 
the onset of saturation corrections as well. Secondly, the argument is based on the principle 
that the mean dipole density will be a good gauge of the saturation corrections. However in 



the case of unitarisation corrections, it has been shown [14] that the use of mean interactions 
can be very misleading. In this section we will perform a more sophisticated analysis of the 
importance of saturation corrections, relative to unitarity corrections. We will also present 
work on understanding wavefunction evolution once saturation becomes a large effect. 



3.1 Measures of saturation 

It is first necessary to develop a measure of saturation. Two methods will be used here. 
Firstly one can probe a wave function with an (unevolved) onium. The two gluon exchange 
amplitude between the wave function and the onium indicates the local density dipoles of size 
similar to the onium and might be similar to the interaction that a dipole in the wave function 
would have with its neighbours. Since the dipole-dipole interaction has spikes wherever two 
sources overlap, one should average the interaction over a region of roughly the onium size. 
The choice of the size of the region over which to average is one of the main uncertainties 
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Figure 5: The probability of a given onium configuration lying in a certain I max , O max bin, 
for y = 10. 

with this method. It also has the disadvantage of being relatively slow, because one needs 
to determine dipole-dipole interactions at many points to determine the saturation over the 
whole extent of the wave function and if one wants the saturation associated with dipoles of 
different sizes, one must use probes of different sizes. 

The second method is more arbitrary, but much faster to compute. Each dipole has an 
"overlap" associated with it. For each pair of dipoles within the onium wave function which 
are separated by less than a certain distance (which will be a function of the dipole sizes) 
the quantity ct|c</c> is added to the overlap for each of the two dipoles (where the smaller 
dipole has size c< and the larger, size c>). This function is loosely based on the fact that the 
total two-gluon interaction amplitude between a pair of dipoles (moving fast relative to each 
other) has the form a|c</c>(l + lnc>/c<)/2 and that it is dominated by a region where the 
two dipoles are close in impact parameter. This method has the advantage that one obtains 
a value of the saturation for each dipole. 

Comparing the two methods one finds that there is a good (linear) correlation between 
them, as can be seen in figure |5| which shows the frequency distribution for configurations with 
various values of I max (the maximum interaction with the probe) and O max (the maximum 
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overlap). One can carry out other comparisons, in particular, considering the location in 
impact parameter of the maximum saturation, and one again finds that the two methods give 
very similar results. In addition, if one examines other values of rapidity, one finds that the 
factor relating I max and O max stays approximately the same. 



3.2 Distribution of I max 

The degree of saturation is obviously just related to the density of dipoles. Therefore the 
distribution of I ma x should be related to the probability of obtaining large densities, which 
has been shown to be exponential [14], the slope of the exponential being proportional to the 
mean density of dipoles. Instead of searching the whole relevant region of impact parameter 
for the location of I ma x > one finds the location of the dipole corresponding to O max and bases 
the search for l max around that point. The probe size is always the same as the onium size. 

One finds, as expected, that I ma x also shows an exponential behaviour, with a slope 
varying in proportion to the local density of dipoles. One also finds that the maximum 
saturation is usually close to the original onium which is also unsurprising since that is where 
the evolution has the longest rapidity range to produce a large number of dipoles. 

While these properties are of general interest, one of the main reason for studying satu- 
ration is to see whether it is likely to affect the onium-onium cross section. So one needs to 
know the typical values of I max and also the correlation between I max for a pair of onium 
configurations and the interaction cross section for that pair of configurations. 

The most appropriate value of rapidity to examine is Y = 12, since this is the rapidity 
for which the mean 1-pomeron amplitude at r = 0, F^\r = 0), is equal to 1 (the unitarity 
bound). Consider first the 'unweighted' curve of figure ||: this corresponds to the probability 
distribution of I ma x in onium-onium collisions with a total rapidity Y = 12. The rapidity 
for each onium is y — — 6. Imax is now defined to be the larger of the Imax 

values 

for the two onium configurations. One sees that typical values of saturation are 0.1 — > 
0.3. The exponential probability distribution for large saturation is relatively clear. More 
relevant though are the distributions weighted with the onium-onium cross section for each 
configuration pair: the weighted distributions correspond to the fraction of the cross section 
coming from configuration pairs with a certain I max , and one sees that the typical saturation 
has increased to be in the range 0.2 — * 0.45. There are two main reasons why configuration 
pairs with larger cross sections are associated with larger saturation. Firstly if one of the onia 
has a high density of dipoles in some region, then the interaction of that region with a second 
onium will be enhanced. A second possibility, which we will consider in more detail later, is 
that configurations which lead to a large interaction (e.g. those with large dipoles) are most 
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Figure 6: The probability distribution of the larg er I max from a pair of onia ('unweighted' 
curve), and the fractions of cross sections coming from configuration pairs with particular 
values of I max ('1-pomeron' and 'unitarised' curves); Y = 12. 



likely to be produced through branching sequences which pass through a high density stage. 

There remains the problem of assessing the effect that this degree of saturation would have. 
One can assume that saturation will start to become important when I max ~ 1 ; but this is a 
very approximate condition, being uncertain to within at least a factor of two. This factor of 
two is critical: if the condition for saturation is Imax ~ 2 then there will be no problem with 
saturation at Y = 12, on the other hand if it is Imax ~ 1/2 then saturation could significantly 
alter the results at this rapidity. If saturation is important then one wants to know how it 
will affect the evolution: for example the saturation might be large in one region of impact 
parameter, but the bulk of the cross section might come from different regions, not affected 
by saturation. The next sections address these and other issues. 

3.3 Saturation in the toy model 

One approach to the problem is again from the toy model, using the condition that the 
scattering amplitude must be independent of the frame in which calculation is performed. 
Equivalently, it shouldn't matter which of the two onia one chooses to evolve. Starting with 
a state containing n dipoles, the rate of branching into a state containing n + 1 dipoles is 
defined to be R n , while the interaction between two states with n and m dipoles respectively 
is F(n,m). If the evolution is partitioned such that for an increase 5Y in total rapidity, the 
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onium with n dipoles evolves by £,5Y and the one with m dipoles evolves by (1 — £)SY, then 
the evolution of the amplitude F will be 



dF 

— = iR n [F(n + 1, m) - F(n, m)\ + (1 - £)R m [F(n, m + 1) - F(n, m)\ (13) 

In the case of one pomeron exchange, one has R n = (ap — l)n and F(n, m) = fa 2 s nm (where 
/ sets the strength of the interaction), giving 

-jy- = K + (1 " 0K«7> - l)fo? s nm (14) 

which is clearly independent of £. 

If for F(n, m) one uses the form which includes unitarity corrections, F(n, m) = [1 — 
exp(— fagnm)], one obtains the following expression for the evolution of the interaction: 



dF s 



£R S n (l - e" /Q S™) + (1 - £)R s m (l - e -/4™)] e -o? s mn ( 15 ) 



dF 

where the superscript s indicates that saturation is being taken into account. The condition 
that the growth of the cross section be independent of £, as well as one's knowledge of the 
branching rate for small n (where saturation should be unimportant) leads to the following 
expression for the rate of branching: 

K = ^^(l-e-f^). (16) 
J a s 

Note that though the interaction with a second onium was used in the determination of this 
branching rate, the result is independent of the state of the second onium. 

To see what effect saturation will have on the final probability distribution of dipoles, one 
can use the following arguments. First take the case without saturation, in the continuum 
limit for n. The probability distribution satisfies 

dP dR n P dnP 

d^ = -^r = - ( ^- 1) ^r- (17) 

Changing variables to Q = nP, ( = {a-p — l)y and v = Inn, one has 

dQ _ dQ 

dc " ~ d^' (18) 

whose solution (from our knowledge that P = exp[— {a-p — l)y — nexp(— (ap — l)y)]) is 

Q = exp - C - e"-<) . (19) 
In the case with saturation one has 
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dy 



dR s n P s 
dn 



(a-p - 1) d(l - e ~ fa s n )P s 



fa s dn 

By an appropriate change of variables, Q s = (1 — e~f a s n )P s /(/a|), ( = (a-p 
v s = ln[(e^ a s" — l)/(/cr|)], one can obtain an equation identical to eq. (|T8|) 



dQ s 



dQ & 



(20) 
l)y and 

(21) 



dC dz^ s 

Because, in the limit of small n and £, Q s ~ Q and i/ s ~ z/, the equation for Q s will have the 
same initial conditions as that for Q and therefore their solutions will be the same. Translating 
this back to the variables P s and n, one has 



P s (n,y) = exp 



-(ap-l)j/ 



(22) 



For e^" 13 " 1 ^, n <C l/(/a|) this approximates to the solution without saturation. The main 
feature of this solution though is that at large y, instead of there being an exponential dis- 
tribution in n with mean exp((a-p — l)y), one now has a solution which has a maximum 
at 



1 



n 



l)y + ln/a| 



(23) 



and which dies off quickly beyond that.0 

One effect of this modified probability distribution is to increase the value of the S'-matrix 
at very large rapidities. If one calculates the S'-matrix in a centre of mass frame (rapidity Y 



divided equally between the two onia) neglecting saturation one obtains [11] 



S{Y) ~ 



In 




/a|e(^- 1 ) y / 2 



(24) 



Note that this is much larger than the value of the S'-matrix that one obtains from typical 
configurations (namely exp(— fa 2 s e( av ~ l ^ Y )). This difference arises because the S-matrix is 
dominated by a small fraction of configurations which have few dipoles. 

Calculating the S-matrix including saturation (done in the lab frame, to reduce sensitivity 
to the use of the continuum limit, though the result is of course independent of frame), one 
obtains 



3 Note that beyond n = {a-p — l)y/fa% the solution eq. (|22|) breaks down because the continuum limit used 
in its derivation is no longer valid (P 3 varies too rapidly over a range Sn ~ 1). Numerical solution shows 
though that there is still a rapid fall off beyond the maximum, having roughly a Poisson distribution, since 
the branching rate n — * (n + 1) becomes independent of n. 



17 



When /a|e^- x ) y / 2 < 1, eqs. (||) and fl2|) give the same result, consistent with the idea that 
saturation should not be important in the centre of mass frame for that energy range. However 
for fale^v-W 2 > 1, eqs. (§|) and (§|) do differ (the saturated case has a logarithmic 
dependence on energy, absent from the calculation without saturation), because saturation 
then matters in the evolution of the onium wave functions, even in the centre of mass frame. 

While this approach works well for the toy model, in the real case, an analogous method 
will turn out to be too complex to be useful. Nevertheless one can also make use of the effect 
of a change of frame on the average (unitarised) amplitude. The simplest pair of frames to 
examine would be centre of mass and lab frames for onium-onium scattering. However one 
finds that, except for logarithmic factors, the 5-matrix behaves in a very similar way in the 
two. Consider instead collisions between a toy onium and a toy "nucleus" consisting of N 
onia. Two frames will be used: a lab frame Lo where the onium is at rest, and a lab frame 
Ltv where the nucleus is at rest. At high energies (neglecting saturation), in frame Lpj, with 
the onium being evolved, the matrix is approximately 

S " = N falser ■ (26) 
Except for logarithmic terms (and the factor N) this has the same leading energy behaviour 
as the centre of mass onium-onium S'-matrix, and is similarly dominated by those rare config- 
urations of the onium which have few dipoles (the common configurations have such a large 
number of dipoles that their contribution to the S-matrix can be neglected). 

In the Lo frame, it is the nucleus which it is evolved. Assuming for the sake of the 
argument that the constituent onia evolve independently (consistent with our neglect of sat- 
uration), configurations with few dipoles will be much rarer, because they require that none 
of the N constituent onia should have evolved much. Accordingly, when one calculates the 
S-matrix, one finds 

So sr j (27) 

which is much smaller than the result in the Ln frame. One way of restoring the corre- 
spondence is by requiring that saturation alter the probability distribution of the number of 
dipoles so that the value of the S'-matrix for a typical configuration of the nucleus is now 

~ <28) 



18 



Equivalently the typical number of dipoles must be of the order of [{a-p — l)Y+m N fa 2 s ]/(fa 2 s ). 
Note that this is just the result of eq. (|23| ) (except for the In N term which is a result of the use 
of nuclei). The advantage of examining the effect of frame changes on the average 5"-matrix 
is that one can transfer the arguments fairly straightforwardly to the case of scattering with 
transverse dimensions. 



3.4 Saturation in the case with transverse dimensions 

The real case, including transverse dimensions, is very much more complicated, to a large 
extent because the dipole picture breaks down at high dipole densities and one needs to take 



into account quadrupoles, hexapoles, and other more complicated colour structures [12, 20]. 
This section will deal with various approaches that one can nevertheless try in order to deduce 
more information about the effects of saturation in the real case, still using the requirement 
of frame independence. 

It is simple to write down the equation which must be satisfied for frame independence 
to hold. This will be the equivalent (with transverse dimensions) of the independence of 
eqs. ( |j~3|) and ( |l5| ) on £. Let 7 be a particular gluon configuration and {7b} be the set of all 
colour configurations with an extra gluon added at transverse position b, with 7b being any 
particular one of these. In the dipole approximation (large N c ), a particular 7b corresponds to 
the gluon having originated from a particular dipole. Define also F~ y (r) to be the interaction 
between configurations 7 and 7' (moving in opposite directions) when they have a relative 
separation in impact parameter of r. Then for the scattering to be frame independent (or 
equivalently for it not to matter which onium is evolved) the rate of evolution ^(7 —* 7b) 
must satisfy 



E I d 2 bi?( 7 - 7b) {F lb>1 , - F 7 ,y) = E / d2hR ^' "> O { F iX ~ F 7,7') (29) 

{7b K {7 b } 

In the case of of the one-pomeron approximation for the scattering amplitude F, and the 
normal dipole evolution kernel for R, if one can show that eq. ( p9| ) holds for a single dipole 
in each onium, then because dipoles in each onium interact (and evolve) independently of 
one another eq. ( |29] ) automatically holds for any collection of dipoles in each onium. The 
fact that eq. (^) holds for the 1-pomeron case is a result of the conformal invariance of 
the dipole-dipole interaction and of the dipole evolution kernel. The details are given in the 
appendix. 

The effect of including multi-pomeron interactions is to connect the interaction of dif- 
ferent dipoles within the same onium, and as a result dipoles no longer necessarily evolve 
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independently of one another: this is the breakdown of the dipole picture of evolution which 
necessitates that one take into account quadrupoles and more complicated terms. Unfortu- 
nately it has not so far been possible to solve eq. ( |29| ) to obtain detailed information on the 
effects of saturation. 

One can nevertheless use frame independence to deduce some limited information. We 
will first examine the effect of saturation on central impact parameters. The argument will be 
along similar lines to the toy-model approach discussed at the end of the last section, except 
that now it will just be sufficient to consider simple onium-onium scattering in lab (L) and 
centre of mass (CM) frames. It has been shown in [ fL4[ that there is a close correspondence 
between toy model and real results for the probability distribution of the number of dipoles 
when this number is high compared to the mean. However the arguments used here depend 
on the probabilities of having unusually low dipole densities, where the toy model is a poor 
approximation to the real case. In a lab frame one finds that the S-matrix is approximately 
S L oc e - c ( y - y °) 2 instead of S L oc e"^"^ in the toy model. Qualitatively, the difference can 
be understood as coming from the extra dynamics associated with the diffusion in dipole sizes. 
In the toy model, the probability of having no evolution is just the exponentially suppressed 
probability of a single dipole not branching over a rapidity Y. In the real case the probability 
(P p (Y)) of a dipole of size b not producing any new dipoles larger than a small size p is 

/ 2asN c Y b\ 
Pp{Y) = exp ^ In -J . (30) 

The extra element which arises is in the choice of p, since in the rapidity interval Y one must 
also require that none of the dipoles produced below size p branch back up into ones of size 
b. Forcing the probability of producing an extra large dipole from the small dipoles to be 
small (approximately e^ ar ~ 1 ^ Y p 2 /b 2 <C 1), one obtains the condition Inb/p oc Y which leads 
to the observed dependence of the S-matrix on Y 2 (since the S'-matrix is dominated by the 
configurations with small interactions and therefore few dipoles). 

The behaviour in the centre of mass frame will be similar, except that the condition for a 
small interaction will be that both onia should be unevolved, so 



S C m oc exp 2 



-c 



2 



exp 



A2 



(31) 



The value of Yq is different from Yq because it depends on the interaction between the two 
onia, the details of which are different in the lab and centre of mass frames. Figure ^ shows 
results from the Monte Carlo simulation, plotting the derivative of InS. The clear straight 
line behaviour at the higher values of Y is the signal of the Y 2 term in the exponent. Fits 
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Figure 7: The logarithmic derivative of the zero impact parameter S-matrix for centre of 
mass and lab frames. Curves under the heading 0-0 are for onium-onium collisions, while 
the 0-N curves are for onium-nucleus collisions, where the nucleus consists of two onia with 
a random relative transverse position of the order of the onium size. 



to the slopes of the derivatives show exactly the predicted factor of two between the lab and 
centre of mass frames. For as — 0.18, the value of c is c ~ 0.042. If one performs similar 
calculations for onium-nucleus collisions, one expects (and sees) very similar behaviour: in 
the Ln frame, the slope is identical to the onium-onium lab frame results, and the Lq slope 
is twice the slope (where the nucleus consists of two onia). 

Since, for onium-onium collisions, the asymptotic value of the S'-matrix is much smaller 
in the lab frame than in the centre of mass frame, one concludes that saturation in the lab 
frame will increase the S'-matrix, and so that the typical value of the S-matrix, including 
saturation will be 

gtyptcal (y) _ e " C ^^ , (32) 

very different from the typical value without saturation, ~ exp[— a s e^ av ~ 1 ^ Y ]. This behaviour 
of the S-matrix can related to the typical fields in an evolved onium or nucleus and this is 
done in |g0| . 
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Figure 8: The 1-pomeron and unitarised scattering amplitudes for onium-onium collisions in 
centre of mass and laboratory frames. 



22 



3.5 Phenomenological modification of the dipole branching 

A second approach to saturation in the real case is intended to give some information on what 
happens at non-zero impact parameters. First one can compare the unitarised amplitude for 
the lab and centre of mass frames as a function of r (the amplitudes have been multiplied 
by r 2 to indicate which regions contribute most to the total cross section). This is shownQ 
in figure ||. The slight mismatch between the 1-pomeron curves at moderate r is due to 
statistical error (the amplitude at larger values of r comes from very rare configurations). As 
one expects, due to the neglect of saturation, the unitarity corrections in the lab frame are 
significantly smaller than in the centre of mass frame. Note in particular though that whereas 
at small impact parameters the difference between the two frames is only a small fraction of 
the unitarisation corrections, at large impact parameters, most of the unitarity corrections 
present in the centre of mass frame are absent in the lab frame. This suggests that in the 
lab frame, saturation effects will in fact be considerably larger than the unitarity corrections 
at large r. One way of investigating this further is to try to include some simulation of the 
effects of saturation. 

The problem is that one doesn't know how to include saturation effects. Our approach will 
be to use a very simple model of saturation which retains the basic properties of normal dipole 
evolution and where saturation changes only the rate of dipole branching, not the spatial 
distribution. This is not expected to give accurate quantitative information on evolution with 
saturation, only to produce some understanding of the types of effects which are liable to be 
important if one does take saturation into account properly. In analogy with the toy model, 
we will use a rate of branching, Rf, of a dipole i of size b into two, of sizes c< (the smaller) 
and c> (the larger) including saturation, of 

Rl{b^ c< , c> ) _f 1 _ e -o i/ n\ n (33) 



Ri(b —> c<,c>) V I ' Oi 

and where Ri is the normal rate of branching (i.e. without saturation). The overlap measure 

of saturation, Oj, is used because it is much faster to calculate, and during the Monte Carlo 

evolution, it will be the calculation of the saturations for each dipole which will be the limiting 

factor. Since the functional form of this saturation is taken from the toy model, one does 

not know what normalisation 0, to use for the overlaps. The procedure used for choosing £1 

will be to find a value which gives a frame independent result for F(r = 0) (r is the relative 

4 Note that the calculation of unitarisation corrections in the lab frame is potentially incorrect, because the 
formula, jf 1 ""') — 1 _ exp(— /), used for the unitarised interaction from a pair of configurations which give 
an interaction / with one pomeron exchange, breaks down when the number of interacting dipoles in either of 
the onia is small and also when the fields in either onium become too strong. 
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Figure 9: The onium-onium scattering amplitude, showing the frame invariance that results 
when both unitarity and saturation corrections are included. 
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Figure 10: Onium-nucleus collisions (nucleus consists of two onia) in various frames: the 
upper set of three curves contains only unitarity corrections, while the lower set has both 
unitarity and saturation corrections. 

impact parameter between the colliding onia) at a single value of rapidity (which is always 
possible), and then to see whether this also works at other values of r and rapidity. This 
has been done for Y = 8, where we take ~ 0.33. Reexpressing this in terms of the a more 
natural quantity, the probe interaction /, this corresponds to a modification of the branching 
rate by a factor ~ (1 + J/1.3), which seems reasonable. 

Figure [9| shows the scattering amplitudes for centre of mass and lab frames, including 
saturation for Y = 12. The one-pomeron and unitarised centre of mass results are also 
plotted for reference. The procedure has worked insofar as the results in the two frames agree 
for small r for the second rapidity. At larger r there is still a small discrepancy between the 
two frames: the lab frame result is too low — there is a little bit too much saturation at 
large r — this discrepancy is however small compared to the overall effect of saturation at 
those values of r. The same is found at all other values of rapidity (including Y = 8). Other 
implementations of saturation (i.e. with different forms for eq. (p3|)), give very similar results, 
as long as the effect of saturation on branching rates is linear for small overlaps. 

As a second check on this form of saturation, one can examine onium-nucleus collisions. 
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As well as the two lab frames discussed earlier we examine also a "centre of mass" frame 
(or more accurately, a frame in which the mean number of central dipoles is the same for 
the onium and the nucleus). Neglecting saturation one expects the Lo frame to have the 
lowest unitarity corrections (since with the nucleus containing all the rapidity, it will have 
the largest saturation corrections), followed by the Ljy and CM frames. Figure |l0| shows this 
ordering of the amplitudes. Once saturation is included (using exactly the same procedure 
as for onium-onium collisions), the three frames give very similar amplitudes, though again 
there is a small discrepancy at larger r. 

It would be interesting to examine how this implementation of saturation affects the typical 
S'-matrix for small impact parameters, and whether this is in accord with the arguments 
presented earlier. However the asymptotic e~ c ( Y ~ Y °^ 2 / 2 behaviour sets in only for relatively 
large rapidities (one needs Y > 14), and calculation of the saturation corrections at these 
rapidities is extremely time consuming (because the calculation time is proportional to the 
square of the number of dipoles) and therefore not feasible. 

Considering the difference between the amplitude in the centre of mass frame with and 
without saturation, at small r, saturation seems to have a fairly small effect on the centre of 
mass amplitude. This is as expected from the basic argument that each onium has half the 
rapidity, and should be well away from saturation when unitarity corrections start to set in. 
At large r, however, the effect of saturation is bigger, and in fact of the same order as the 
unitarisation correction (which relates to the observation for figure || that at larger r there is 
a big difference between the centre of mass and lab frame results, indicating that saturation 
must have a significant effect at large r). 

The dynamics involved in producing the effects saturation at large r are complex: because 
the unitarisation corrections are small, the densities of dipoles must be small, and one might 
be misled into thinking that saturation corrections should unimportant. One reason why 
this argument fails is that there are a variety of ways in which large r interaction can occur. 
A large dipole can be produced early in the evolution, which can then evolve to produce 
many other dipoles of similar size and the final density of dipoles will be indicative of the 
relevant saturation corrections (these evolution paths tend to also have significant unitarity 
corrections). Another evolution path involves only small dipoles for most of the rapidity 
and a branching to larger dipoles at a late stage. The evolution of the small dipoles may 
have been subjected to significant saturation; but once the branching occurs to large dipoles, 
unitarity corrections and any further saturation corrections (at the large scale) are no longer 
affected by what is happening on the small scale. The history of saturation can therefore have 
an important effect on the evolution, even though the diluteness of the final configuration 
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Figure 11: Saturation at different stages of evolution. The maximum dipole size at a given 
rapidity is c max , and O max is the maximum overlap of any dipole between c max and c max /3 
(essentially a measure of the saturation on the scale of c max ). Errors are large at large r 
because of a lack of statistics in that region. 
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might make one think one could neglect saturation. This is illustrated in the upper part of 
figure [ll]: at two values of rapidity (y = 6 and y = 8), the maximum dipole size, c max (y) is 
determined. The maximum overlap, O max is then calculated for all dipoles with size c, where 
Cmax > c > c max /3, in effect giving the saturation on the scale of the largest dipole. The 
figure shows the ratio O max (y = 8)/O max (y = 6), the increase in saturation (on the largest 
scale), as a function of c max (y = 8)/c max (y = 6), the increase in the largest dipole size. The 
main point is that when there is a large increase in dipole size, the ratio of the saturations 
drops to below one, confirming that if one looks only at the final configuration, the history of 
the saturation can be hidden. A similar reduction of saturation with increasing rapidity can 
probably also take place for the small dipoles produced from large ones (since the small ones 
will be more dilutely spread out). 



A second effect is shown in the lower part of figure 11: the production of a large dipole 
tends to be a rare event. By a simple probabilistic argument, the larger the number of dipoles 
on a small scale, the greater the likelihood that at least one of them will branch to give a 
large dipole, which means that evolution paths that produce a large dipole are more likely to 



include an intermediate high saturation stage. The lower plot of figure 11 demonstrates this 
by showing the maximum saturation (again on the largest scales) at y = 4 as a function of 
the ratio of the largest scales at y = 8 and y = 4. As the ratio of the scales increases so does 
the maximum intermediate saturation. 

Finally in the centre of mass frame, one has the possibility of situations where one of the 
onia is small and moderately dense, while the other has a dilute distribution of small dipoles 
(the only ones which will interact with the small dense onium) at large r. The saturation 
corrections in the evolution of the small dense onium will then be of the same order as the 
unitarity corrections between the two at large r onia, i.e. there is an asymmetry between the 
two onia (and arguments for neglecting saturation when unitarity corrections set in, rely on 
the symmetry between the two onia). This is distinct from the r = case where moderate 
densities of dipoles in each onium led to strongly unitarised amplitude. 

The conclusions therefore, from the simple simulation of saturation effects, and from the 
other arguments presented above, is that at large impact parameters, it is not necessarily safe 
to neglect saturation corrections compared to unitarity corrections, and further that large 
dipoles, even when dilute, may have been produced via a high saturation stage. This might 
be of some concern when determining the corrections to the total cross section, which has 
a significant component coming from moderate to large r. However, to the extent that it is 
possible to gain any quantitative information from our simple simulation of saturation, one 
finds that in the total cross section, saturation corrections (up to Y ~ 16) are no more than 
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a third of the unitarity corrections. 
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Appendix 

This section demonstrates the invariance of the one-pomeron scattering amplitude on the 
choice of frame. We show that it is invariant for a single dipole in each onium. Because each 
dipole evolves independently, and because its interaction with the other onium is independent 
of the other dipoles in the onium to which it belongs, this will imply that the onium-onium 
scattering amplitude is invariant under changes of frames at the one-pomeron level. 

Let the first dipole have sources at b\ and 62, and the second dipole have sources at 63 and 
64 (with the b n being complex variables, and b n their complex conjugates). The condition for 
invariance of the scattering amplitude on the division of a small amount of rapidity between 
them is, eq. (|29|), which can be written as 
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To see that this holds one notes that each of the integrands (including the db^db^ measure) 
is invariant under conformal transformations. The particular conformal transformation 



b + A 
Cb-1 



(35) 



with 



A = 6 1 -6 2 + ft 3 -& 4 (36) 
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6163 - 6 2 6 4 
6I 1 " 6 2 1 + 63 1 - 6J 1 



b 2 v - &r v 

has the property of mapping b\ — > 63, 62 — ► 64, 63 — > b\ and 64 — > 62, thus converting the lower 
integral in eq. ( |35|) to the upper integral (and vice- versa). Since the integrands are invariant 
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under conformal transformations, this proves the equality, and therefore the invariance of 
the one-pomeron onium-onium scattering amplitude on the frame in which the calculation is 
performed, if the evolution is carried out using the dipole kernel and the interaction is by two 
gluon exchange. 
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